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Abstract. In this talk 1 recent developments of the theory of preheating after inflation 
are briefly reviewed. In inflationary cosmology, the particles constituting the Universe 
are created after inflation due to their interaction with moving inflaton field (s) in the 
process of reheating. In inflationary models motivated by supergravity, both bosons 
and fermions are created. In the bosonic sector, the leading channel of particle pro- 
duction is the non-perturbative regime of parametric resonance dominated by those 
bosons which are created exponentially fast with the largest characteristic exponent. 
In the fermionic sector, the leading channel corresponds to the regime of parametric 
excitation of fermions, which respects Pauli blocking but differs significantly from the 
perturbative expectation. In supergravity we also have to consider production of grav- 
itinos and moduli fields, which are cosmologically dangerous relics. We discuss the 
derivation of the gravitino equations in curved space-time with moving background 
scalars. We describe recent results on the production of gravitinos from preheating, 
which may put strong constraints on the inflationary models. 



Preheating after Inflation 

According to the inflationary scenario, the Universe initially expands quasi- 
exponentially in a vacuum-like state without entropy or particles. At the stage 
of inflation, all energy is contained in a classical slowly moving inflaton field (p. 
The fundamental Lagrangian £(<p, x, ip, Ai, h^, ...) contains the inflaton part with 
the potential V(<f>) and other fields which give subdominant contributions to gravity. 
The Friedmann equation for the scale factor a(t) and the Klein-Gordon equation 
for cj)(t) determine the evolution of the background fields. In most of the models, 
soon after the end of inflation, an almost homogeneous inflaton field 4>(t) coherently 
oscillates with a very large amplitude of the order of the Planck mass around the 
minimum of its potential. This scalar field can be considered as a coherent super- 
position of inflatons with zero momenta. The amplitude of oscillations gradually 
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decreases not only because of the expansion of the universe, but also because en- 
ergy is transferred to particles created by the oscillating field. At this stage we shall 
recall the rest of the fundamental Lagrangian which includes all the fields interact- 
ing with inflaton. These interactions lead to the creation of many ultra-relativistic 
particles from the inflaton. Gradually, the inflaton field decays and transfers all 
of its energy non-adiabatically to the created particles. In this scenario all the 
matter constituting the universe is created from this process of reheating. If the 
creation of particles is sufficiently slow, the particles would simultaneously inter- 
act with each other and come to a state of thermal equilibrium at the reheating 
temperature T r . This gradual reheating can be treated with the perturbative the- 
ory of particle creation and thermalization. However, typically particle production 
from coherently oscillating inflatons occurs not in the perturbative regime but in 
the non-perturbative regime of parametric excitation. Indeed, let us consider a 
simple toy model of chaotic inflation with the quadratic potential V((f>) = \m^ 2 
and C int = — \g 2< t> 2 X 2 describing the interaction between the inflatons and other 
massless Bose particles x- The quantum scalar field x in a flat FRW background 
has the eigenfunctions Xk(t) e~ ikx with comoving momentum k. The temporal part 
of the eigenfunction obeys the equation 

Xk + 3- X k+ [%-iR + g 2 <nxk = v (l) 

a \a z J 

with vacuum-like initial conditions: Xk — in the far past. The coupling to the 
curvature £R will not be important in the presence of the interaction (but would 
lead to gravitational preheating in the absence of the interaction). In this model, the 
inflaton field <p{t) coherently oscillates as <p{t) « <&(t) sin (m^t), with the amplitude 
= ^|=- ^ decreasing as the universe expands. The smallness of g 2 alone does 
not necessarily lead to the perturbative excitation of Xk modes. To check whether 
the interaction term g 2 <f) 2 in eq. (1) is perturbative or not, it is convenient to use 
a new time variable z = mt and the essential dimensionless coupling parameter 
1 = Scalar metric fluctuations in this model are compatible with cosmology if 
the inflaton mass is m ~ 10 _6 M P ; therefore, it is expected that q ~ 10 10 g 2 ^> 1 for 
not negligibly small g 2 . In fact, a consistent setting for the problem of x-particle 
creation from the 0-inflaton requires q ^> 1 even without additional assumptions 
about g 2 [2]. Indeed, it is known that if we have two scalars and x, then the 
latest stage of inflation will be driven by the lightest scalar. The square of the 
effective mass of the x-field includes a term g 2 <f) 2 . Inflation is driven by the 0-field 
if its square mass m 2 is smaller than g 2 <f) 2 . This leads to the condition q 1, i.e. 
to the creation of x-particles in the resonance regime. 



Supergravity and the Early Universe 



To make the next step beyond toy models of particles interactions with the 
inflaton, we have to choose the fundamental Lagrangian £(4>, x, i>, Mi Kk, •••)• We 



may expect that the low-energy physics of the early universe will be described by 
the general four-dimensional N — 1 supergravity-Yang-Mills-matter theory [1]. A 
rather lengthy N = 1 phenomenological supergravity Lagrangian 2 begins with the 
terms 



A particular choice of the form of the Lagrangian is motivated and notations are 
given in [3]. In Eq. (2) we choose the minimal Kahler potential K = j^r, where 

is the complex conjugate of $j. The last term of the 1st line is the scalar 
potential V(§i). The equations of motion based on the first line should describe 
inflation, which is a challenging problem by itself. For simplicity we will take a 
prototype model of the superpotential W = \/\& 3 /3, which for |$| <C Mp leads 
to the effective potential A<£> 4 . We will illustrate some effects of preheating in 
supergravity with this model. 



Let us briefly recall the basics of the bosonic preheating. Consider the creation of 
X-particles due to the g 2 x 2( P 2 interaction. To get a feeling for the way in which x~ 
particles are created from inflaton oscillations, we have to understand the character 
of solutions of eq. (1) for the mode functions Xkif). For preheating of bosons we 
need the first line of the Lagrangian (2). Consider the model with the inflaton 
potential V((j>) = \\<t> A in an expanding universe. The problem of particle creation 
in this theory can be reduced to a similar problem in Minkowski space-time. This 
can be realized with the conformal transformation of the scalar field — > tp = acp 
and with the conformal time variable r = \f\(p J -%r. Therefore, \\<j) 4 theory is 
sometimes dubbed the conformal theory. In conformal variables, the Klein-Gordon 
equation for y?(r) is reduced to an equation in flat space-time. Its solution is 
ip(r) ~ (p /(r), where the amplitude of the oscillations <p is constant until the 
backreaction of created particles is taken into account. The time- dependence of 
the oscillations in this theory is not sinusoidal, but given by an elliptic function 
/(r) = an (t, -^j. Eq. (1) for quantum fluctuations Xk can be simplified in this 
theory. Using a conformal transformation of the mode function X^t) = a(t)xk{t) 
in Eq. (1) we obtain 




X 3 W - X 1 TPXi - (e K/2 WW Xt Xj + h.c) 

+ ( V/^Vvn^Vw? + ^Lp&YXi + *P R x t e K/2 vw + h.c. 



) 



+ ... 



(2) 



Preheating of Bosons 



2 ) When I show a viewgraph with the full supergravity Lagrangian at the cosmology conferences, 
somehow people start to laugh hysterically. 



X£+( K 2 + g/ 2 (r))x fe = , 



(3) 



where k 2 = and q — The equation for fluctuations does not depend 
on the expansion of the universe and is completely reduced to the similar prob- 
lem in Minkowski space-time. This is a special feature of the conformal theory 
\\<p A + \g 2 4> 2 X 2 - The mode equation (3) belongs to the class of Lame equations. 
The combination of parameters q = g 2 /X ultimately defines the structure of the 
parametric resonance in this theory. This means that the condition of a broad para- 
metric resonance does not require a large initial amplitude of the inflaton field, as for 
the quadratic potential. The strength of the resonance depends non-monotonically 
on the value of the q parameter. The stability/instability chart of the Lame equa- 
tion (3) in the variables (k 2 , was constructed in [6]. To see how the general 
theory works, let us consider parametric resonance of inflaton fluctuations fluctua- 
tions 4>k due to the self-interaction of inflaton field in A0 4 theory. Using conformal 
transformation of the mode functions y?fc(r) = a(t)(f>k(t), the equation for ipk can be 
reduced to the general equation (3) with the particular value q = 3. The equation 
for fluctuations Xk when q = 3 can be solved analytically [6]. The resonance in the 
"inflaton" direction is weak; the maximal value of the characteristic exponent of 
the fluctuations 0^ oc e MlT is n ~ 0.036. Let us however consider a supersymmetric 
version of the conformal theory. In the first line of the Lagrangian (2) we will put 
the superpotential W = \/A$ 3 /3, which at <C M P gives us a scalar potential with 
the conformal properties V(<&) = |($ 2 ) 2 . In supersymmetric theories, all scalar are 
complex. The scalar field $ will have two components (0,0). Let us assume that 
inflaton direction corresponds to the real component 0, Re $ = A=. and initially 

= 0, Im $ = 0. The equation for the mode function of fluctuations in the direction 
can be obtained with the conformal transformation </?fc(r) = a(t)(pk(t) and can 
be reduced reduces to the general equation (3) with the parameter q — 1. Again, 
in this case the problem can be solved analytically [6]. The factor q — 1 instead 
of 3 in eq. (3) makes a big difference, which manifests the subtlety of the para- 
metric resonance. The resonance in the direction is much stronger and broader 
than the resonance in the inflaton direction, fe oc e^ T , with p, ps 0.147. Thus, the 
supergravity generalization of bosonic preheating in the conformal theory makes a 
big difference in the bosonic preheating due to the self-interaction. Note that the 
character and strength of the parametric resonance for the general equation (3) 
depends on the shape of the effective inflaton potential V{4>). Moreover, the inves- 
tigation of the resonance in an expanding universe typically cannot be reduced to 
the study of the regular stability/instability chart. If theory is not conformal, (say 
due to the mass term m 2 2 in V ((/>)) the parameter q in an expanding universe is 
time-dependent. For the broad resonance case g > 1 this parameter can jump over 
a number of instability bands within a single oscillation of the inflaton field, and the 
concept of stability/instability bands is inapplicable here. Parametric resonance in 
this case is a stochastic process [5]. 



Parametric Excitations of Fermions 



A simple model for the inflaton's interaction with Fermi particles x is a Yukawa 
term hx4>X- F° r instance, consider second line of eq. (2) and ignore mixing between 
X and t/jfj, (which corresponds to the rigid SUSY limit). For our toy model with W = 
!v^V£ 3 and <C M P: the "mass" term of the chiral fermion x is equal to \/2X(j), 
which corresponds to the inflaton-fermion interaction \f2\x4>X- F° r fermions, the 
Pauli exclusion principle prohibits the occupation number from exceeding 1. For 
this reason, it has been silently assumed that fermions are created in the three-legs 
perturbative process <fi — > XX where individual inflatons decay independently into 
pairs of ^-particles. Let us, however, consider the Dirac equation for a massless 
quantum Fermi field x(t,%)'- 

[^V, + h<P(t)]x = 0, (4) 

where V M is the derivative with the spin connection. We are using the representa- 
tion of gamma matrices where 70 = diag(i,i, —i, —i). Here, similar to the bosonic 
case, the inflatons producing fermions also act not as individual particles but as a 
coherently oscillating field <f>(t). Let us consider more general model |A0 4 + hip<j)ip. 
This is a conformal theory in the sense that the problem of fermion production by 
the inflaton in an expanding universe can be reduced to equations in Minkowski 
space-time. Indeed, let us perform a conformal transformation of the involved fields, 
(p = a<f> and \f = a 3 ^ 2 x, and use a conformal time variable, r as in the previous 
section. The equation for the eigenfunctions of the quantum fluctuations in this 
theory can be reduced to a second-order equation for an auxiliary field X(r,x), 
so that \I/ = [7 M V M + hip] X. The eigenmodes of the auxiliary field have the form 
-Yfc( r ) e+ik ^Rri where the R r are eigenvectors of the Dirac matrix 7 with eigen- 
value + 1. The temporal part of the eigenmode obeys an oscillator-like equation 
with a complex frequency which depends periodically on time 

XZ+(K 2 + qf-i^f)x k = 0. (5) 

The dimensional comoving momentum k enters the equation in the combination 
^2 = k 2 ; therefore, the natural units of momentum are y/Xif. The background 
oscillations enter in the form /(r) of the previous section, having unit amplitude. 
The imaginary part of the frequency in Eq. (5) guarantees the Pauli blocking for 
the occupation number n k . The results for n k can be formulated as follows [4]. 
Even though the Yukawa interaction contains a small factor h, one cannot use the 
perturbation expansion in h. This is because the frequency of the background field 
oscillations is proportional to another small parameter y/\. The combination of 
the coupling parameters ^ = q ultimately determines the strength of the effect. 
The growth of fermionic modes occurs in the non-perturbative regime of parametric 
excitation. The modes get fully excited with occupation numbers n k — 1 within tens 
of oscillations of the field 0, and the width of the parametric excitation of fermions 



in momentum space is about q 1 / A \/\(j) . F° r instance, in the case of h = \/2A, 
q = 2, the modes will be excited in about ten oscillations, and the width will be 
about \f\<pQ. 

Equations for the Gravitino 

The rest of my talk is based on the recent project [3] 3 . Let us now consider the 
third line in the supergravity Lagrangian (2), which describes the gravitino field ip^. 
In a general background metric and in the presence of complex scalar fields with 
non-vanishing VEV's, the equation for the gravitino has on the left hand side the 
kinetic part = e^^'-y^Vpip^, and a rather lengthy right hand side. We will use 
the long derivative with the spin connection and Christoffel symbols, for which 
T^nlv = 0. Apart of varying gravitino mass m, the right hand side contains a chiral 
connection and various mixing terms like those in the 3rd line of (2). For a self- 
consistent setting of the problem, the gravitino equation should be supplemented 
by the equations for the fields mixing with gravitino, Xi from (2), as well as by the 
equations determining the gravitational background and the evolution of the scalar 
fields. Let us make some simplifications. We consider the supergravity multiplet 
and a single chiral multiplet containing a complex scalar field $ with a single 
chiral fermion \- This is a simple non-trivial extension which allows us to study 
the gravitino with a non-trivial FRW cosmological metric supported by the scalar 
field. A nice feature of this model is that the chiral fermion \ can be gauged to 
zero so that the mixing between ip^ and \ i n (2) is absent. We also can choose 
the non- vanishing VEV of the scalar field (inflaton) in the real direction, as in the 
previous sections. First we will derive the equation for a spin 3/2 field in a curved 
background metric with non-vanishing VEVs for the scalar fields. 

1 Spin 3/2 Field Equations in External gravitational and Scalar 

Fields 

From (2) we can obtain the equation for the gravitino field (in case of single chiral 
multiplet x an d vanishing Im $) = m( y T)^ fJ,u i/j u , where ^ u = \ (7^7^ — ^ v ^) 
and gravitino mass m = m(0(r)) is given by m = e K l 2 j^r. This equation can be 
transformed into the form 

(in \ 
v f ,-- lll j 1 ^u. (6) 

The gravitino equation (6) is a curved spacetime generalization of the familiar 
gravitino equation (^+mo)Vv = in a flat metric, where m is a constant gravitino 
mass. The generalization of the first constraint equation d^ip^ = can be obtained 
from the equality 7 M -R M = 27 /i "X> /x ^„ and reads 
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V^- ty-f^ + Irn-y"^ = . (7) 

The generalization of the second constraint equation 7 /1 ^ = can be obtained 
from the equality V^W = — \G vp '~) u ip p (dropping the torsion term) and is 

^mVVv + (7^m) 7 ^ - (fymty" = —G^V • (8) 

In our case m = m(r) and from (8) one can find an algebraic relation between 7 O, 0o 
and y-ipi: 

7 >o = h% • (9) 

Here A is a matrix which will play a crucial role in our description of the interaction 
of the gravitino with the varying background fields. If p and p are the background 
energy-density and pressure, we have G® = Mp 2 p, G\ = —Mp 2 pS l k , and one can 
represent the matrix A as follows: 

, p-3m 2 Mp 2m'a~ 1 Mp 

a= ^w p + ^ttw p =Ai+j ° A2 - (10) 

We shall solve the equation (6) using the constraint equations in the form (8) and 
(9). We use a plane- wave ansatz %p p ~ e lkx for the space-dependent part. Then 
ipi can be decomposed into its transverse part ipf, and to the longitudinal part ipf 
which is defined by the trace Y^i- Two degrees of freedom of ip^ are associated 
with the transverse part ipj , which correspond to helicity ±3/2 and two degrees of 
freedom are associated with (or ip ) which correspond to helicity ±1/2. For the 
helicity ±3/2 states we have to derive the equation for ipf. We apply decomposition 
ip = ipj ± ipf to the master equation (6) for p, — i and obtain 

( 7 "^ + Ya 1 " + ma ) ^ = ° ' (U) 

The transformation ipf = a~ 1//2 \I/f reduces the equation for the transverse part to 
the free Dirac equation with a time- varying mass term ma, c.f. eq. (4). In the 
previous section we explained how to treat this type of equation. The essential 
part of tyf is given by the time-dependent part of the eigenmode of the transversal 
component X T (r)), which obeys second-order equation (c.f. (5)): 

X't + (k 2 + (ma) 2 - i{ma)') X T = . (12) 

The corresponding equation for gravitino with helicity 1/2 is more complicated. We 
have to find k l ipi and 7V* ■ The equation for the components k % ipi can be obtained 
from the constraint equation (7). The equation for ^ipj can be derived from (6). 
Using ipg from (9), we get an equation for 7^ 



(d v + B-tk- ll0 A) 7 Vi = 0, 



(13) 



where B = — ^A — ^7 (1 + 3A The time-dependent factor of the spinor 7*^1, 
which we denote as fk{ r )-, obeys a second-order differential equation. By the sub- 
stitution f k (r) = E(r))X L (rj), with E = (-A*) 1 ' 2 exp (- J v drj ReB) , the equation 
for the function fk(r) is reduced to the final oscillator-like equation for the time- 
dependent mode function X L (rj): 

X'i + (\A\ 2 k 2 + Q 2 L - iQ' L ) X L = 0. (14) 

Here a -1 !^ = |9 T ln A* + y~A 2 + \ma (—1 + 3Ai). For an arbitrary background 
FRW metric a(r) and background scalar field 0(r), equation (14) may lead to ill- 
defined physics. For instance, if \A\ > 1, it describes noncasual propagation of spin 
3/2 particles. In the context of preheating, the background fields 0(r) and a(r) are 
oscillating. Naively, one could expect that \A\ is oscillating too. In this case one 
would reach a pathological conclusion that the strongest parametric excitation of 
spin 3/2 eigenmodes will be with the highest momenta k. Noncausality and other 
defects of spin 3/2 fields interacting with an external electromagnetic field are well 
known [7]. It is interesting to compare the equations for a spin 3/2 field in an 
external electromagnetic field and in external gravitational/ scalar fields. Then we 
will show how the apparent defects of (14) are resolved in supergravity. 

2 Spin 3/2 Field Equations in an External Electromagnetic Field 

The equations for a charged spin 3/2 field interacting with an external electro- 
magnetic field A^ with field strength tensor in flat space-time were derived in 
[7]. After simple manipulations, they can be re- written in the form of the equation 
of motion 

fir^n + m Vv = (A* _ ^T 7 ^) 7 ^ ' ( 15 ) 
and two constraint equations 

2>^_ + ^ mo7 ^ = , (16) 

^7^ = _Ii>^, (17) 

where = e / -f 5 e flupcr Fp a . Here T>^ = <9 M — ieA^. The system of equations (15), 
(16), (17) for flat case with EM background is similar to the equations (6), (7), (8) 
for cosmological problem up to definition of derivatives instead of V^, a constant 
mass m instead of the altering mass m(r), and F pu instead of G pv . It is interesting 
that a spin 3/2 field interacting with external fields can be described in the unified 
way. In [7] it was shown that for some configurations of the electromagnetic field, 
the propagation of spin 3/2 field violates causality. 



3 Self- consistent Gravitino Problem in Supergravity 



Contrary to the inconsistent setting of the problem of a spin 3/2 field in an 
arbitrary EM or gravitational field, the spin 3/2 gravitino field in supergravity 
should be consistent. In this section we will show that the gravitino equations in 
an expanding universe with moving scalars are consistent. We will use a model 
with a single chiral multiplet. Let us concentrate on the matrix A given by (10). 
In the models where the energy-momentum tensor is determined by the energy 
of a classical scalar field and $ depends only on time we have p = |$| 2 + V, 
p = |$| 2 - V. The scalar potential is V($) = e K \VW\ 2 - 3m 2 Mp. Also, we have 
m! = ae K / 2 UW&/Mp. Therefore, the matrix A can be rewritten in terms of $ and 
e K l 2 VW only 

|<l>| 2 - e K/2 VW 2 2$e K / 2 VW 

A = rrro , ^7772 + 7o ,-, , , ^ TT „ 2 • (18) 



|$| 2 ± \e K ' 2 VW\ 2 |$| 2 ± \e K ' 2 VWf 

From this form of A it follows that 

\A\ 2 = A\ + A 2 2 = l (19) 

for an arbitrary superpotential W. Thus A can be represented as A = 
— exp (2i Jl^ dt n{rj) ). Using the Einstein equations, one obtains p = WW + A, 

where the correction A = 0(Mp 1 ) is given in [3]. The expression for fi becomes 
much simpler and its interpretation is more transparent if the amplitude of oscil- 
lations of the field $ is much smaller than Mp. In the limit $/Mp — > one has 
fj, = d<$,d<$,W . This coincides with the mass of both fields of the chiral multiplet (the 
scalar field and spin 1/2 fermion) in rigid supersymmetry. When supersymmetry 
is spontaneously broken, the chiral fermion, goldstino, is 'eaten' by the gravitino 
which becomes massive and acquires helicity ±1/2 states in addition to the helicity 
±3/2 states of the massless gravitino. With this form of A the gravitino equation 
(14) becomes 

X'[ + {k 2 + Q 2 L - iVL' L ) X L = . (20) 

with a~ 1 fli — A* — |i?sin2/ \idt — \m{\ ± 3 cos 2 / fxdt). Equation (20) is consis- 
tent and describes the creation of gravitinos from preheating. The solution of the 
consistency problem in the equation (20) is that the symmetries imprinted to the 
Supergravity provide \A\ = 1. But it does not mean that A is a constant. The 
matrix A does not become constant even in the limit Mp — > oo. The phase of A 
rotates when the background scalar field oscillates. The amplitude and sign of A 
change two times within each oscillation. Consequently, the relation between 7 o, 0o 
and Y^Pi a l so oscillates during the field oscillations. This means that the gravitino 
with helicity 1/2 (which is related to ^o) remains coupled to the changing back- 
ground even in the limit Mp — > oo. In a sense, the gravitino with helicity 1/2 



remembers its goldstino nature. This is the main reason why gravitino production 
in this background in general is not suppressed by the gravitational coupling. The 
main dynamical quantity which is responsible for the gravitino production in this 
scenario will not be the small changing gravitino mass m(t), but the mass of the 
chiral multiplet /J,, which is much larger than m. 

Problem of Gravitino Over-Production from Preheating 

As an example, consider the model with the superpotential W = \/A<I> 3 /3. The 
parameter \x for this model is given by /i — V2X(f). It rapidly changes in the interval 
between and v / 2A0o- Initially it is of the same order as H and to, but then H and 
m rapidly decrease as compared to /i, and therefore the oscillations of \x remain the 
main source of the gravitino production. In this case the production of gravitinos 
with helicity 1/2 is much more efficient than that of helicity 3/2. The theory of 
production of gravitinos with helicity 1/2 in this model is similar to the theory of 
production of spin 1/2 fermions with mass \/2A0 by the coherently oscillating scalar 
field in the theory A0 4 /4, which we considered above. Growth of helicity 1/2 grav- 
itino modes (20) occurs in the non-perturbative regime of parametric excitation. 
The modes get fully excited with occupation numbers — 1 within about ten 
oscillations of the field 0, and the width of the parametric excitation of fermions in 
momentum space is about \/A0o. This result violates the cosmological constraints 
on the abundance of gravitinos with mass ~ 10 2 GeV by 4 orders of magnitude 
[3]. The most dangerous gravitino over-production (by 14 orders of magnitude) 
occurs in the class of inflationary models where V(<fi) does not have a minimum 
and where preheating is gravitational (NO models) [2]. Thus the investigation of 
the non-thermal gravitino production in the early universe may serve as a useful 
tool helping us to discriminate among various versions of the cosmological theory. 
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